Coulomb-enhanced resonance transmission of quantum SINIS junctions 
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Coherent charge transfer through a baUistic gated SINIS junction is mediated by the resonant 
tunneling via the Andreev states. Extra charge accommodated on the Andreev levels partially com- 
pensates the charge induced by the gate voltage preserving the electron wavelength and maintaining 
the resonance conditions in a broad range of gate voltages. As a result, the transparency of the 
junction as well as the supercurrent trough it can be substantially increased as compared to the 
zero-Coulomb case. 

PACS numbers: 74.45.-l-c, 74.50.-|-r,73.63.-b 



According to the classical picture of the Coulomb 
blockade in gated structures [l[ the Coulomb interaction 
usually suppresses electronic transport through small 
conductors by introducing an additional energy barrier 
associated with the charging energy. In this Letter we 
show that in ballistic SINIS junctions (here S stands for 
a superconductor, N is a normal-metal island, and I is 
an insulator), the Coulomb interaction can rather stimu- 
late the supercurrent than suppress it. The supercurrent 
through a junction containing a ballistic normal island is 
mediated by Andreev states with energies controlled by 
the effective transparency of the double-barrier structure 
which is determined by the wave length of electrons. The 
extra charge on the dot affects not only the energy but 
also the wave length of electrons. This makes the trans- 
mission extremely sensitive both to the gate voltage and 
to the charge on the Andreev states localized in the nor- 
mal island. We show that charging of Andreev states 
can actually preserve the condition of resonance tunnel- 
ing: the charge adjusts itself so as to compensate the 
deviation of the chemical potential caused by the change 
in the gate voltage thus increasing the transparency of 
the device for the supercurrent. 

The .system. - We consider the normal island in 
the form of a short single-mode ballistic conductor 
(wire) connected to bulk superconducting leads via low- 
transparency contacts. This corresponds to the common 
experimental situation which is realized, in particular, in 
recent experiments on nanowires 0]. The length d of a 
conductor is much less than the superconducting coher- 
ence length, ^, i.e., |A| <C 5E where 5E ^ hvp/d is the 
level spacing in the double-barrier structure, vp is the 
velocity of electrons in a quantum conductor, and |A| is 
the superconducting gap in the leads. The short con- 
ductor (sometimes called an Andreev quantum dot) is 
weakly coupled to the superconducting leads and forms 
a low-proximity device. This setup is similar to that con- 
sidered in Refs. Q where the charging of Andreev levels 



was studied in the limit of small Coulomb interaction. 

To characterize the efficiency of the Coulomb interac- 
tion in the SINIS we introduce the dimensionless param- 
eter 7 = Ec/SE ~ e^d/ChvF, where Ec = e^/2C is 
the Coulomb energy of the island, C being its capaci- 
tance. The maximum value of 7 is less than e^/hvpx, 
with the island capacitance C ~ xd\n.(d/a), where a is 
the transverse dimension of the wire and >c is the effective 
dielectric constant in the presence of a substrate. For a 
typical value vp ~ 10* cm/s the corresponding 7 is small; 
thus the assumption 7 ^ 1 is adequate for practical de- 
vices. The same parameter characterizes the Luttingcr- 
liquid effects which thus can be neglected. However, 
in spite of the small value of 7, the Coulomb energy can 
still be either larger or smaller than the superconduct- 
ing gap. To account for nonlinear back-action effects 
on the contact transparency associated with charging of 
the Andreev level in case of a considerable Coulomb en- 
ergy Ec ~ |A|, we adopt a mean- field model similar to 
that employed in Ref. |5[ for normal junctions. In this 
respect, our model differs from those considering trans- 
port through quantum dots having a single fixed electron 
level 0. 

Model. - The charging Hamiltonian has the form 
H,^Ec{N-No)\ N^y^l ii{r)i;^{v)d\, (1) 

the sum being taken over the spin indices a and the in- 
tegration domain is the volume V of the island. Here 
iVo is the background charge supplied by ions and by the 
external circuit. The kinetic energy has the usual form. 
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where ^(r) = I[6{x + d/2) + S{x - d/2)] describes the 
insulating barrier between the island and the leads, x is 
the coordinate along the conductor. 
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Superconductivity in the leads is described by the 
usual BCS Hamiltonian where A(r) = |A|e*-*^^'^ inside 
the leads, i.e., for x < —d/2 ot x > d/2, respectively, 
while A = in the normal conductor, —d/2 < x < d/2. 

In this model, the states of the system are spin- 
independent. The advanced and retarded Green func- 
tions in the real-frequency representation can be found 
from the Dyson equation 



ydVL(ri,r)Gf(-4)(r,r2) = (5(ri-r2) 



(2) 



where 

L(ri, r) = Go\ri)S{r - n) + 2Echiri)K{ri,r) , 

are matrixes in the Nambu space, and 

Ho = -n^vVSm -Ef + V{r) + Uch{v) , 



Uc = 2Ec{N - No), N =- d-'rTvKir, r) . 

Jv 

The function h{r) = 1 if r belongs to the island and 
h{r) =0 otherwise. The kernel 



(3) 



is determined by the Keldysh function, Gf^ (ri,r2) ~ 
[Gf (ri,r2) - G^^(ri,r2)] /i(e). Here /i(e) = 1 - 2n, 
and is the occupation number. In equilibrium rig is 
the Fermi function, and /i(e) = tanh(e/2T). The energy 
Uc describes the variation of the charge on the island 
due to the change in the number of excitations while N 
has a meaning of the excitation charge on the island. 

Bogoliubov-de Gennes (BdG) equations. ~ Let us ex- 
pand the retarded and advanced Green functions in a set 
of orthogonal and normalized functions as 

= iVn)- The functions u and v satisfy the 
linear equations with effective potentials determined by 
Kuu{r,ri) and if™(r,ri): 



HoUn{r) +A{r)vn{r) + Ech{r) /^^^(r, ri)u„(ri)dVi 

Jv 

-Ech{r) Kuv{r,ri)vn(ri)d^ri ^ EnUn{r) , (4) 
Jv 



Hovnir) -A*(r)w„(r) +Ech{v) A':„(r, ri)w„(ri)d^ri 

Jv 

+Ech{r) /v*„(r,ri)u„(ri)(iVi (r). (5) 

Jv 



The kernels i^'„„(r, ri) and ii'„„(r,ri) are determined 
self-consistently. They can also be expanded into the 
same wave functions. 
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The kernels obey the symmetry if„ij(r, ri) = 
-K*yir,ri) = A"*„(ri,r) and /f™(r,ri) A"*„(r,ri) = 
Kuv , I") due to the symmetry of the functions 
ue — > , ve ^'^*-E with respect to £" — > —E. The 
excitation charge is also expressed through m„ and w„, 
^=iE„^n: where 



dV K{r)u„ir) - <(r)z;„(r)] h{E,,) (7) 



is the number of excitations on the state n, it determines 
the average excitation charge Qn = eNn carried by the 
state. The sum runs over the states with both positive 
and negative energies. It can be written as a sum over 
only positive-energy states: TV = J2n e„>o ^n- 

Andreev states. - The normal conductor has a single 
transport mode parameterized by the coordinate x along 
it. The eigenstates in a short low-transparency junction 
have the form of narrow particle-like and hole-like reso- 
nances such that only two of them fit into the sub-gap 
range of |A|. Due to proximity to the superconducting 
leads these states transform into the Andreev states lo- 
calized over the distances of the order of ^ near the junc- 
tion. All other resonances are essentially the same as in 
the normal state; thus they can be ignored if one is in- 
terested in energies of the order of the superconducting 
gap. The two Andreev states will be labeled with n = +1 
and n = — 1. One of them has a positive and other has 
a negative energy Ei ~ —E^i > 0. In what follows we 
assume that the transparency is so small that T <C d/£^. 
In this limit the Andreev states are mostly concentrated 
inside the normal conductor. We will also consider ener- 
gies much smaller that the superconducting gap because 
it is this energy range where the supercurrent enhance- 
ment is most pronounced. In this case the kernels Eqs. 
([6|) in Eqs. ^ and ([5|) contain contributions only from 
the two Andreev states, n = ±1. Due to the symmetry 
with respect to i? — )- —E, 

Kuuix^xi) ^ [ui{x)ul{xi) - vl{x)vi{xi)]fi{Ei) , 
Kuvix,xi) ^ [ui{x)vl{xi) + vl{x)ui{xi)]fi{Ei) . 



Denoting ui 



u e and vi = i;+e*P"^ 



V e and neglecting the difference between the mo- 
menta of u and v in the kernels, pu = Pv = P (for a short 
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junction, pu — Pv ^ E/vp ^ l/<^) wc find 

l.d/2 



-d/2 

d/2 
-d/2 



dxi [Ki,uix,xi)ui{xi) - Kuy{x,xi)vi{xi)] 

= -u^{x)N^-vl{x)Ah, 
dxi [Kl^{x,xi)vi{xi) + Kl^{x,xi)ui{xi)] 
= -vi{x)Ni +ul{x)Mi , 



Ni being the number of excitations on the state n = 1 
while Ml = 2 j'^'l^^^ui{xi)vi{xi) fi{Ei) dxi. 

The BdG equations (jlj, ([5|) inside the normal conduc- 
tor take the form 



where the effective chemical potential is 

fi = Ef -Uc + EcNi . 
For u,v* (X e*P^' we find from here 



p'/2m^^i^±jEf-E^„ 



(8) 



(9) 



where Eg — Ec\Mi\ is the Coulomb gap. Using the 
solutions of the BdG equations without the Coulomb in- 
teraction, one can show that Mi ~ T. The Andreev 
state energy is (see below) Ei > {hvF/d)T- Therefore, 
the Coulomb gap Eg ^ EcT in Eq. ([9]) can be neglected 
when 7 <C 1. In this limit Eq. ([9]) reduces to its usual 
form, p^ /2m — ji ~ ±-Ei, while the functions u and v 
inside the normal conductor are decoupled. 

The BdG equations for the entire system can thus be 
written in a standard form 



Ho-h{r)EcNi ui{r) + A{r)vi{r)=EiUi{r), (10) 



Ho - h{r)EcN, 



wi(r)+A*(r)ui(r) =£;iz;i(r). (11) 



The BdG equations (|TOl) . (|TT|) determine the energy of 
the Andreev bound states as a function of the phase dif- 
ference (j) = xr ~ XL- For low transparency, T ^ d/^, 
the energies E+i = ±{hvp/ d)eA in the range £' ^ |A| 
have the form IS, M 



tA = V^2+7-2cos2(0/2) 



(12) 



where T ^ Ti/2 is the transparency of two contacts in 
series. The quantity z measures the shift of the chemical 
potential /i inside the conductor, Eq. ([5]), from one of 
the resonances. It is defined through the equation 

PFd/h + r] + {d/hvF)[EcNi - Uc] ^ TTk + z (13) 



where r] = arctany^ is the scattering phase shift and 
k is an integer. The range of applicability of Eq. (|T2|) is 
EA < d/i. 

As it follows from Eq. 0, C/c = C/c + '2'EcNi where 
U'(j = 2Ec{N' - No) and N' = J^n^i ^n- Let us denote 



y ^ PFd/h + T] — dU'fjIhvF — T^k . 

The quantity y depends on the gate voltage through 
the background charge A^o- The excitation number A^' 
changes little if the chemical potential varies within one 
resonance, z < I. Therefore, the chemical potential shift 
y can be considered as a parameter controlled by the gate 
voltage. 

For a localized state (sl 



fdEi 



h{Ei). 



(14) 



The subscript (j) means that the phase is kept constant. 
As it follows from Eq. ([8]) for the effective chemical po- 
tential of the Andreev state, dfi/dz = hvF /d. Therefore, 
Eq. ((HI) yields 



y = z + ^^^fi{Ei) 
oz 



(15) 



We consider zero temperatures and put fi{Ei) = 1. 
This can be done if temperature satisfies T ^ T6E ~ 
AT{£,/d). Equation ([T5|) is the central result of this pa- 
per. Equations (|12p and (|15p determine eA as a function 
of y and Eq. then determines the charge A^i on the 
Andreev level. 

Results. - Equations (fT^ and ([T5|) yield 



(l + 7/eA)v'4-r2cos2(0/2) 



(16) 



We see that at 7 7^ ( i.e., due to the Coulomb inter- 
action) there appears a substantial region of gate volt- 
ages, y 7, where energy is close to its smallest pos- 
sible value, Tcos((/)/2). This is more pronounced if the 
Coulomb energy is larger than the gap Ec/\A\ ^ 1 (es- 
timates will be given later). Indeed, since ea ^ d/^ the 
ratio 7/eyi > £'c/|A|. Yet more drastic effect is expected 
if T ^ 7. In this case the term 7/eA in brackets in Eq. 
([T6| dominates. As a result, for y < 7 



r cos((/)/2) 



(17) 



According to Eqs. (fT6|) or (fTT]). ca and z remain small 
eA,z<^^ in the range of gate voltages y < 7. Thus the 
energy satisfies Ei <C Eq, |A|, and the Kondo regime is 
not realized. 

Figure[l] (upper panels) shows the solutions of Eqs. (fT6|) 
for T = 0.1 and for two values of the Coulomb interac- 
tion, 7 = 0.2 (left panel) and 7 = (right panel). One 
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observes that e^i — J- as approaches tt if y < 7. The 
decrease in energy is, of course, cut off when |(/) — 7r| ^ 
(r/|A|)(d/^T) because of the decreasing factor in 
Eq. (jlSp . This energy behavior has a crucial impact on 
the supercurrent given by the expression Q : 



2e fdEi 



h{Ei 



Jo- 



Tsin ( 



(18) 



where jo = evpT/'^d. Since the energy remains small in 
a wide range of gate voltages, the supercurrent is essen- 
tially enhanced as compared to its zero- Coulomb value 
for the same gate voltage. In particular, for compara- 
tively large Coulomb interaction, 7 2> T, we have from 

Eq. inD 

J /jo = 2v/l-yV7' sin(0/2) . 



The critical current jc — 2jo\/l — ll^ is reached at 
= TT. It is roughly by factor 7/T larger than its value 
without the Coulomb interaction, jc = joT /y, reached 
at 7r/2 for the same gate voltage y ^ T ■ The su- 
percurrent Eq. (fTS]) is shown in Fig. [1] (lower panels) for 
T = 0.1 and for two value of the Coulomb interaction, 
7 = 0.2 (left panel) and 7 = (right panel). One ob- 
serves significant enhancement of the supercurrent for a 
given gate voltage. 

Discussion. - Wc have shown that charging of the 
Andreev level leads to a strong enhancement of the super- 
current in a wide range of gate voltages (parameterized 
by the quantity y) . The charge of the Andreev state com- 
pensates the deviation of the chemical potential caused 
by the change in the gate voltage thus preserving the res- 
onance transmission of the double barrier structure and 
increasing the transparency of the device. 

The Coulomb interaction is characterized by the di- 
mensionless parameter 7 ~ dEc/hvp- The upper esti- 
mate gives 7 ~ /hvp ^ 1- Our approximation requires 
small 7. Nevertheless, the factor 7/T that determines the 
enhancement of the supercurrent can reach quite large 
values due to low transparency. A substantial enhance- 
ment of supercurrent occurs already when Eq ^ |A|. A 
typical gap A (as for Al) is about 1 K. The Coulomb en- 
ergy Ec = e^/2C of 1 K corresponds to the capacitance 
C « 8 • 10~^ cm. Therefore, to have Eq ^ A the char- 
acteristic size of the normal conductor should be smaller 
than ~ lO"'^ cm which can be easily realized in practice. 
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FIG. 1. Andreev state energy (upper panels) and the super- 
current (lower panels) in the presence of Coulomb interaction, 
7 = 0.2 (left panels) and without it, 7 = (right panels) as 
functions of the phase difference <j) and the gate voltage y for 
T = 0.1. The curves for 7 = and 7 = coincide at y — >■ 0. 
The curves for the supercurrent at 7 = 0.2 for y = and 
y — 0.01 overlap, while at 7 = they significantly difi^er at cf> 

close to TT. 
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